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We give a classification of spherically symmetric kinematic self-similar solutions. This 
classification is complementary to that given in a previous work by the present authors 
[Prog. Theor. Phys. 108, 819 (2002)]. Dust solutions of the second, zeroth and infinite kinds, 
perfect-fluid solutions and vacuum solutions of the first kind are treated. The kinematic self- 
similarity vector is either parallel or orthogonal to the fluid flow in the perfect-fluid and 
vacuum cases, while the 'tilted' case, i.e., neither parallel nor orthogonal case, is also treated 
in the dust case. In the parallel case, there are no dust solutions of the second (except when 
the self-similarity index a is 3/2), zeroth and infinite kinds, and in the orthogonal case, there 
are no dust solutions of the second and infinite kinds. Except in these cases, the governing 
equations can be integrated to give exact solutions. It is found that the dust solutions in 
the tilted case belong to a subclass of the Lemaitre-Tolman-Bondi family of solutions for the 
marginally bound case. The flat Friedmann-Robertson- Walker (FRW) solution is the only 
dust solution of the second kind with a — 3/2 in the tilted and parallel cases and of the 
zeroth kind in the orthogonal case. The flat, open and closed FRW solutions with p — — fl/3, 
where p and /i are the pressure and energy density, respectively, are the only perfect-fluid 
first-kind self-similar solutions in the parallel case, while a new exact solution with p = fj,, 
which we call the "singular stiff-fluid solution", is the only such solution in the orthogonal 
case. The Minkowski solution is the only vacuum first-kind self-similar solution both in the 
parallel and orthogonal cases. Some important corrections and complements to the authors' 
previous work are also presented. 

§1. Introduction 



Scale-invariance is one of the most fundamental properties of classical gravity. 
In the vacuum case without a cosmological constant, the gravitational constant G 
is the only dimensional constant that appears in the field equations in Newtonian 
gravity, while there appears another physical constant, the speed of light, c, in general 
relativity. Neither a characteristic time nor a length scale can be constructed from 
them. This implies that there exists a class of scale-invariant solutions, which are 
often referred to as self-similar solutions. There also exist self-similar solutions in 
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scalar-tensor theories of gravity. (See Ref. 1) for an investigation of the perfect-fiuid 
case in Brans-Dicke theory.) In contrast, there may not exist self-similar solutions 
in quantum gravity, since, by definition, there appears a third dimensional constant, 
the Planck constant, h, with dimensions independent of that of G and c. Matter 
fields that include some dimensional constants could also lead to violation of the 
scale-invariance. In general relativity, matter fields consistent with self-similarity 
include no constants with dimensions independent of those of G and c. 

Self-similar solutions have been studied extensively, especially in spherically sym- 
metric systems. This is not only because the self-similarity ansatz leads to the re- 
duction of the governing equations from partially differential equations to ordinary 
differential equations, whose mathematical treatment is generally simpler, but also 
because there exist many applications. Although self-similar solutions are only spe- 
cial solutions of the field equations, they often play an important role as an attrac- 
tor or an intermediate attractor. In particular, there is a self-similarity hypothesis 
which states that solutions in a variety of astrophysical and cosmological situations 
may naturally evolve into self-similar forms even if they are initially more compli- 
cated. 2 ) In recent years, the validity of this hypothesis in spherically symmetric 
gravitational collapse has been reported both in Newtonian gravity 3 ) " 9 ) and general 
relativity. 10) " 15) 

Self-similarity in general relativity is defined by a self-similarity generator, which 
is a vector field. It includes three subcases, called the "parallel," "orthogonal" and 
"tilted" cases. They are distinguished by the relation between the generator and a 
timelike vector field, which is identified as the fluid flow, if it exists. The "tilted" 
case, in which the generator is neither parallel nor orthogonal to the fluid flow, is 
the most general case among them. 

Self-similarity in general relativity was first defined by Cahill and Taub as the ex- 
istence of a generator called a homothetic Killing vector field. 16 ) Such self-similarity 
is often called homothety. For the perfect-fluid case, homothety restricts a barotropic 
equation of state to take the linear form p = K/j,, where p and \i are the pressure 
and energy density, respectively, and K is a constant. The dominant energy con- 
dition requires — 1 < K < 1. Spherically symmetric homothetic solutions in the 
tilted case have been intensively investigated in various situations. It has been 
found both in Newtonian gravity and general relativity that a homothetic solution is 
an attractor solution in gravitational collapse 3 ) " 15 ) or a critical solution associated 
with critical phenomena. 17 )~ 20 )' 14 )' 21 ) Spherically symmetric homothetic solutions 
have been used in modeling the growth of primordial black holes 22 )' 23 ) and the 
evolution of a cosmic void. 24 ) Those in the tilted case were classified for the dust 
(K = 0) case by Carr 25 ) and for the perfect-fluid case by Carr and Coley. 26 ) (See 
also Refs. 27) -29).) Coley showed that the Friedmann-Robertson- Walker (FRW) 
solution is the only spherically symmetric homothetic perfect-fluid solution in the 
parallel case. 30 ) Mcintosh showed that a stiff fluid (K = 1) is the only compatible 
perfect fluid with homothety in the orthogonal case. 31 ) 

One possible type of self-similarity compatible with more general equations of 
state is kinematic self-similarity. It is one natural generalization of homothety and 
is defined by a generator called a kinematic self-similarity vector field. 32 )' 33 ) Kine- 
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matic self-similarity is classified with respect to an index into four kinds, the first, 
second, zeroth and infinite kinds. In the context of kinematic self-similarity, homo- 
thety is considered the first kind. Thus, when we also take the relation between a 
self-similarity generator and a fluid flow into consideration, there are twelve subcases 
of kinematic self-similarity. Although there are some works on kinematic self-similar 
solutions, 34 ) " 37 ) the importance of kinematic self-similarity is still less clear than 
that of homothetic solutions. In a previous work by the present authors 39 ) (I), it 
is assumed that the equation of state is polytropic or p = K/j, and a classifica- 
tion of spherically symmetric kinematic self-similar solutions was presented. 38 )' 39 ) 
In that work, we only treated the cases of a non-zero pressure and vacuum, leav- 
ing out the pressureless (i.e. the dust) case. Spherically symmetric dust solutions 
can be obtained in closed form with two arbitrary functions and are known as the 
Lemaitre-Tolman-Bondi family of solutions. 40 ) Among them, spherically symmetric 
homothetic dust solutions have been used in analyzing naked singularities, 41 ) " 43 ) in 
modeling the growth of primordial black holes 22 ) and in describing the evolution of 
a cosmic void. 24 ) Spherically symmetric kinematic self-similar dust solutions could 
have as many important applications as homothetic solutions. In this paper, we 
extend I, add some corrections and complements, and complete the classification. 

The organization of this paper is the following. In §2, the basic equations in 
a spherically symmetric spacetime are presented, and kinematic self-similarity is 
briefly reviewed. In §3, we classify the spherically symmetric kinematic self-similar 
dust solutions of the second, zeroth and infinite kinds. In §4, we give corrections and 
complements to 1. 39 ) In §5, we classify the spherically symmetric homothetic perfect- 
fluid solutions together with vacuum solutions for which a homothetic Killing vector 
field is parallel or orthogonal to the fluid flow. Section 6 is devoted to summary and 
discussion. The flatness of 3-surfaces for kinematic self-similar dust solutions are 
presented in Appendix A. We adopt units for which c = 1. 

§2. Spherically symmetric spacetime and kinematic self-similarity 

The line element in a spherically symmetric spacetime is given by 



where df2 2 = dO 2 + sin 2 Odip 2 . We consider a perfect fluid as a matter field, 



where is the four- velocity of the fluid element. Here we have adopted comoving 
coordinates. Then, the Einstein equations and the equations of motion for the perfect 
fluid are reduced to the following simple form: 



(2-1) 



= P(t, r)g^ u + [v(t, r) + p(t, r)]U^U u , 
U^ = (-e*, 0,0,0), 



(2-2) 
(2-3) 



4 



H. Maeda, T. Harada, H. Iguchi and N. Okuyama 



m r = Air[xR r R , 
mt = —4irpR t R 2 , 

= -R tr + $ r Rt+&tRr, 
1 

2G 

where the subscripts t and r indicate differentiation with respect to t and r, and 
m(t, r) is the so-called Misner-Sharp mass. The auxiliary equations from the Einstein 
equations are as follows: 



-Ml + e- 2 *R t 2 



e- 2 *R r 2 ) 



(2-6) 
(2-7) 
(2-8) 

(2-9) 



'W 



R 



+ 2&i 
R 



+ e 



2<P-2'P 



2%-2%<Zv + 
R R 



R r 
~R 



-87rG/ue 

,2^ 



2<P 



-SiiGpe 2 *, 



Rtt R+ — j- 
i? R 



Rt\ 
RJ 



R r 
~R 



2 Rr T 



(2-10) 



(2-11) 



Five of the above eight equations are independent. 

On a hypersurface of constant t, the two-sphere R = 2Gm is a marginally 
trapped surface (or an apparent horizon). The region R < 2Gm is trapped, while 
the region R > 2Gm is not trapped. If we assume R r > 0, outgoing (ingoing) null 
geodesies satisfy 

^ = ±e*-*, (2-12) 
dr ± y ' 

where the subscripts + and — denote outgoing and ingoing, respectively. Along the 
trajectory of an apparent horizon, the relation 

2 



/ds 


) =-e 2 * 


7 dt 




\dr 


AH/ 


\dr 


ah) 



— e 



2&-2<P 



= —e 



2<P 



AH 



dt 
dr 



(2-13) 
(2-14) 



dt dt 
dr AH dr 

is satisfied, where the subscript AH indicates that the derivative is taken at an 
apparent horizon. It is found from Eq. (2T4) that the trajectory of the apparent 
horizon is spacelike when 

dt 
dr 



dt 
dr 



dt 

AH < dr 



(2-15) 



is satisfied on it, while it is timelike when 



dt 
dr 



dt 

AH < ~dr 



or 



dt 
dr 



dt 

+ dr 



AH 



(2-16) 



is satisfied on it. The trajectory of an apparent horizon is an outgoing (ingoing) null 
when 



dt 
dr 



AH 



dt 
dr 



(2-17) 
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is satisfied on it. 

A kinematic self-similarity vector £ satisfies the conditions 

£^ = aU„, (2-19) 

where /i^ = g^v + UfjU v is the projection tensor, denotes Lie differentiation along 
£, and a and 5 are constants. 32 )~ 34 ) The similarity transformation is characterized 
by the scale- independent ratio a/S, which is referred to as the similarity index. 

In the case 5 ^ 0, 5 can be set to unity, and the kinematic self-similarity vector 
£ can be written 

^-^ + «I +r l- (2 - 20) 

if it is tilted, i.e. neither parallel nor orthogonal to the fluid flow. In the case a = 1, 
which corresponds to self-similarity of the first kind (for which (3 can be set to zero), 
it follows that £ is a homothetic vector, and the self-similarity variable £ is given by 
r/t . In the case a = 0, which corresponds to self-similarity of the zeroth kind (for 
which (3 can be rescaled to unity), the self-similarity variable is given by 

£ = re' 1 . (2-21) 

In the case a / 0, 1, which corresponds to self-similarity of the second kind (for 
which (3 can be set to zero), the self-similarity variable is given by 

<=(5^- (2 - 22) 

In the case (5/0, self-similarity implies that the metric functions can be written 

R = rS(0, $ = ¥ = (2-23) 

In the special case that 5 = and the self-similarity is referred to as that of 

the infinite kind (here, a = 1 is possible). The kinematic self-similarity vector £ can 
be written 

In this case, the self-similarity variable is given by 

£ = ~ t - (2-25) 

In the case (5 = 0, self-similarity implies that the metric functions can be written 

R = S(0, e # = e^), e* = e*K>/r. (2-26) 

It is noted that £ is a Killing vector in the case that 5 = and a = 0. We omit the 
bars on # and \P in §3.1 for simplicity. 
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If the kinematic self-similarity vector £ is parallel to the fluid flow, it can be 
written 



e dx» ~ f{t) dt' 



(2-27) 



where f(t) is an arbitrary function, and the self-similarity variable is now r. When 
we set f(t) = t, the metric can be written 

ds 2 = _ t 2( a -l) e 2*(r) dt 2 + t 2 dr 2 + ^(^2^2 (2 . 2g) 

for the first, second and zeroth kinds, and 

ds 2 = _ e 2*(r) dt 2 + dr 2 + S ( r f dQ 2 ( 2 -29) 

for the infinite kind. We omit the bars on 3 in §3.2 and §5.2 for simplicity. 

If the kinematic self-similarity vector £ is orthogonal to the fluid flow, we have 



(2-30) 



where g(r) is an arbitrary function, and the self-similarity variable is then t. When 
we set g(r) = r, the metric can be written 



„2a 



dt 2 + e 2 ^dr 2 + r 2 S(t) 2 df2 2 



for the first, second and zeroth kinds, and 

e 2*(t) 



ds 2 = -r 2 dt 2 + 



dr 2 + S{tfdQ 2 



(2-31) 



(2-32) 



for the infinite kind. We omit the bars on & in §3.3 and §5.3 for simplicity. 

§3. Spherically symmetric kinematic self-similar dust solutions 

Setting Pi = P2 = in the basic equations for a perfect fluid given in 1, 39 ) we 
obtain those for a dust fluid. Since we classified the vacuum case in 1, 39 ) we omit 
the vacuum case in this section. 

3.1. Tilted case 

In the tilted case for dust solutions of the second, zeroth and infinite kinds, 
= is obtained from basic equations. Therefore we set exp(^) = Co, where Co is 
a positive constant, in this subsection. 

3.1.1. Self-similarity of the second kind 

In the case of self-similarity of the second kind, the Einstein equations imply 
that the quantities m and fx must be of the forms 

J2 



2Gm = r 



8irGfi = 



(3-1) 
(3-2) 
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where £ = r/(at) 1 / a . The basic equations are the following: 



1V1\ -+- iW-^ 


w cfir Q I c'^ 

— Wio l^o ■+■ <5 J, 


(o-o 


1 A/f„ _L fl/f' 
oiv±2 ~r 


— VV 2 D [J -r <-> ), 


(1. A 




n 


(O-O 




n 

— u, 


/'Q £ 
(O-O 




— Q\~\ „— 2^/ c i q'\2i 

— £[l — e (5 + i J J, 


[6-1 


a M 2 


C c'2 — 2 


(3-8 


w[s 


= -W^'S + 2S'), 


(3-9 


{2aW 2 + W 2 \)S 


= -W 2 (V'S + 2S'), 


(3-10 


S" + S' 


= (S + S')&, 


(3-11 


S'(S' + 2&S) 


= a 2 W 2 S 2 c 2 , 


(3-12 


2S(S" + 2S') - 2&'S(S + S') 


= -S' 2 -S 2 + e 2 *(l-W 1 S 2 ), 


(3-13 


2S(S" + aS') + S' 2 


= 0, 


(3-14 


{S + S'f 


= e 2 ^, 


(3-15 



where the prime denotes differentiation with respect to ln£. Here we have implicitly 
assumed at > 0. When we define £ = r/(— at) 1 /" for at < 0, basic equations 
identical to Eqs. (3-3)-(3T5) are obtained. We extend solutions analytically to the 
region in which at < if they are analytic at t = 0. Such extended solutions can be 
obtained when we allow solutions to be defined for at < 0. 

Equations (3-15) and (3-7) give M\ = 0, which implies W\ = from Eq. (3-3), 
because \S + S'\ = exp(^) cannot be zero. Equation (3-14) gives the evolution 
equation for S, 

2y' + 3y 2 + 2ay = 0, (3-16) 
where y = S' / S. This equation can be integrated to yield 

S=\ Sl -s C a \ 2/3 , (3-17) 

where sq and si are integration constants. In the present case, sq / must be 
satisfied, because so = gives a vacuum spacetime. The relation M 2 = SS' 2 /(a 2 c^) 
is obtained from Eq. (3-8), while W 2 = (3 — 2a)y 2 /[a 2 c^(\ + y)\ is obtained from 
Eqs. (3-12), (3-11) and (3-14). The resulting solution is 

p 2<p - r 2 

e — Cq, 

c 2* _ [fi + ((2Q/3) - l)s r a ] 2 
\si - so£~ a | 2/3 
[sir Q + ((2a/3) - l)s at] 2 

r 4a/3| Sir a _ SQa t\2/3 

R = r\si- s C a \ 2/3 

= r 1 ~ 2a / 3 \ Sl r a -s at\ 2 / 3 , 



(3-18) 

(3-19) 
(3-20) 
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4(3 - 2a)s§f 



-2a 



9c§t 2 (si - s C- a )[si + ((2a/3) - l)s r a ] 

4(3 - 2a)s 2 l a 2 

9c^(sir a - s at)[sir° + ((2a/3) - l)s atf] ' 



(3-21) 



2Gm - ^H— r 2a 



^W *~2a 



(3-22) 



When the time-reversal solution is also considered, we can choose so > and s\ > 
without loss of generality. With the coordinate transformations t — > i = cot and r — > 

obtain a one-parameter family of solutions parametrized 
by re = (sxCq/ (s a))\c / (s a)\ 2a ^ 3 ~ 2a ^ for each a, 

_ 9(^/(3-^) + (2a/3 - l)t) 2 

+ (3-2a)2| K f3«/(3-2a) _t|2/3 ^ 

+ f 2| Kf 3a/(3-2a)_^|4/3 dfi 2 j (3 . 23) 

~ r „ = 4(3 - 2a) 

M 9( K f 3Q /( 3 - 2 ") + (2a/3-l)F)(Kf 3Q /( 3 - 2Q ) -t)' 1 ' 

2Gm = jjf 3 . (3-25) 
9 

This solution was obtained previously by Carter and Henriksen 32 ) and has been 
treated by several authors. 35 )' 37 ) This solution belongs to the Lemaitre-Tolman- 
Bondi family of solutions for the marginally bound case. We call this solution the 
second-kind kinematic self-similar Lemaitre-Tolman-Bondi (KSS-LTB) solution. We 
omit the bars on i and f hereafter for simplicity. If s± = (i.e., re = 0), this solution 
gives the flat FRW solution for any a / 0,1. When s± / (i.e., re / 0), a ^ 3/2 
must be satisfied for the solution to be a non-vacuum solution. This solution with 
re is asymptotically flat FRW solution for r — ► oo with a<0or3/2<a and for 
r — > with < a < 3/2 (a / 1). It is noted that for re / 0, r — > +0 corresponds to 
R — > +oo for a > 3/2 and to i? — > +0 otherwise. The energy density is positive in 
the regions satisfying 

t < Kr 3a /( 3 ~ 2a \ [3/(3 - 2a)]rer 3Q /( 3 ~ 2a ) <t for < a < 3/2, a + 1 (3-26) 

and 

t < [3/(3 - 2a)]rer 3a /( 3 - 2a ), Kr 3 «/( 3 - 2 ") <t for a < 0, 3/2 < a. (3-27) 

The singular surfaces are t = Kr 3a K 3 - 2a ) and t = [3/(3 - 2a)]rer 3a /( 3 ~ 2a ) . The 
former is a shell-focusing singularity, because the circumferential radius R is zero 
there, while the latter is a shell-crossing singularity, because the quantity 

f {(3/(3 - 2a))rer 3a /( 3 ^ 2a ) - t} / ( K r 3a / ( 3 ~ 2 ^ - t) 1 ' 3 
dR_\ for t < Kr 3a /( 3 ~ 2a \ 

dr " {t - (3/(3 - 2a))rer 3a /( 3 " 2 «)}/(t - « r 3a/(3-2a))i/3 l* 28 ) 

I for t > Kr 3a /( 3 - 2a ) 
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is zero there. For the trajectory of the shell-crossing singularity, the relation 



R = 



2a 



3 -2a 



2/3 



3/(3-2q) 



(3-29) 



is satisfied. The fluid is collapsing in the region satisfying t < nr 3a ^ 3 2a ) and 
expanding in the region satisfying t > Kr 3Q /( 3 - 2a ) ) because the quantity 

^ = -sign(Kr 3a /( 3 - 2Q ) - t )\r\Kr^/^- 2 ^ - t\^ 3 (3-30) 

is negative in the former case and positive in the latter case. The trapped region is 
given by 



KT 



3o/(3-2q) 



(8/27)r 3 < t < Kr 3«/(3-2a) 



and 



Kr 3a/(3-2a) < f < Kr 3a/(3-2a) + ( 8/27 ) r 3. 

The apparent horizons are 

t = K r 3a/{3 ~ 2a) + (8/27)r 3 , (AH1) 

and 

t = Kr 3«/(3-2a) _ ( 8 / 27 ) r 3. (AH2) 



(3-31) 
(3-32) 
(3-33) 
(3-34) 



Here we consider the trajectories of the apparent horizons. We consider only the 
regions with jjl > in which R r > are satisfied. From Eqs (3-33) and (3-34), we 
obtain 



dt 
dr 



__ \ {3qk/(3 - 2a)}r( 5Q - 3 )/( 3 " 2a ) + (8/9)r 2 for AH1, 
AH ~ 1 {3qk/(3 - 2a)}r( 5Q - 3 )/( 3 - 2a ) - (8/9)r 2 for AH2. 



(3-35) 



AH1 has its extremum at r = n = [27aK/{8(2a - 3)}] (3 ~ 2q)/{9(1_q)} for 3/2 < a 
or a < 0, while does not for < a < 3/2 (a ^ 1). AH2 has its extremum at 
r = r 2 = [27ok/{8(3 - 2q)}]( 3 - 2 «)/{9( 1 -«)} for < a < 3/2 (a / 1), while does not 
for 3/2 < a or a < 0. For the null geodesies on the apparent horizon, the relation 



dt 
dr 



± 



' ±|{3qk/(3 - 2a)}r( M )/( 3 ~ 2a ) - (4/9)r 2 
±|{3qk/(3 - 2a)}r( 5a - 3 V( 3 - 2a ) + (4/9)r 2 | 



for AH1, 



(3-36) 



for AH2 



are satisfied. First we consider r > 0. For a < 0, 3/2 < a, the trajectory of AH1 is 
spacelike (timelike) for 



r 9(l-a)/(3-2a) < (>) _ 2 7q K /[2(3 - 2a)], 



(3-37) 
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while it is outgoing null for 

r 9(l-a)/(3-2a) = _ 27aK /[2(3 - 2a)] EE ^(l-<*)/(3-2a) . (g.gg) 

For < a < 3/2 (a / 1), the trajectory of AH1 is timelike. For a < 0, 3/2 < a, the 
trajectory of AH2 is timelike. For < a < 3/2 and a^l, the trajectory of AH2 is 
spacelike (timelike) for 

r 9(l-a)/(3-2a) < (>) 2 7q K /[2(3 - 2a)], (3-39) 
while it is ingoing null for 

r 9(i-«)/(3-2a) = 27a K /[2(3 - 2a)] = r »(i-°)/(3-2«) (3 . 40) 
Next we consider r = 0. For both AH1 and AH2, the relation 

^ = ±sign(a/(3 - 2a)) (341) 

is satisfied for a < 1 and 3/2 < a, so that both AH1 and AH2 are ingoing (outgoing) 
null at r = for a < and 3/2 < a (0 < a < 1). For 1 < a < 3/2, the left-hand-side 
of equation (3-41) is equal to ±2 and =f2 for AH1 and AH2, respectively, so that 
they are not null at r = 0, and the results for r > can be applied for r = 0. 

Schematic figures for this solution are displayed in Figs. 1-4, along with those 
for its time-reversal solutions. For all a (7^ 0, 1 or 3/2), the solution is unphysical 
in region II, in which the energy density is negative. We consider only regions with 
positive energy density. 

For a > 3/2, region I on the left in Fig. 1 represents a decaying white hole, while 
region I in its time-reversal solution represents a growing black hole. Region III on 
the left in Fig. 1 represents the outer region of a contracting shell-crossing singularity 
that starts from i? = +ooatt = — 00. The singularity reaches the physical center 
at t = — 0. Region III in its time-reversal solution represents the outer region of an 
expanding shell-crossing singularity that starts from R = +0 at t = +0. 

For < a < 3/2 and a 7^ 1, region I on the left in each of Figs. 2 and 3 
represents the inner region of an expanding shell-crossing singularity, which forms 
at the physical center at t = +0. Region I in each of the time-reversal solutions in 
Figs. 2 and 3 represents the inner region of a contracting shell-crossing singularity 
that starts from i? = +ooatt = — 00. The singularity reaches the physical center 
at t = —0. Region III on the left in each of Figs. 2 and 3 represents the formation 
of a shell-focusing singularity from a regular spacetime. Because the singularity is 
covered for 1 < a < 3/2, this solution represents the formation and growth of a black 
hole. For < a < 1, further study is needed to determine whether this singularity 
is naked or covered. Region III on the right in each of Figs. 2 and 3 represents a 
decaying white hole. The white hole singularity disappears at t = and leaves a 
regular spacetime. 

For a < 0, region I on the left in Fig. 4 represents a decaying white hole, while 
region I in its time-reversal solution represents a growing black hole. Region III on 
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Region I (Rr>0, U >0) 

Region II (R<0, U <0:unphysical) 

Region III (RrX), U >0) 




Fig. 1. The figure on the left is for the second-kind KSS-LTB solution [Eqs. (3-23)-(3-25)] with 
k 7^ for a > 3/2, and that on the right is for its time time-reversal solution. AH, SFS 
and SCS indicate an apparent horizon, shell-focusing singularity and shell-crossing singularity, 
respectively. Each diagram is divided by SFS and SCS into three regions, which we call regions 
I, II and III. Region II is unphysical, because fi < there. A solid curve represents an apparent 
horizon whose trajectory is timelike, while a dashed curve represents an apparent horizon whose 
trajectory is spacelike, r = r\ corresponds to the extremum of AHl. The trajectory of AHl is 
outgoing (ingoing) null at r = r$ in the left (right) figure. It is noted that r — > +0 corresponds 
to R — > +oo for this a. 

the left in Fig. 4 represents the formation of a non-central shell-crossing singularity 
from a regular spacetime. The singularity forms at R = +oo at t = +0 and contracts 
with time. For t > 0, it represents the inner region of the singularity. Region III 
on the right in Fig. 4 represents the inner region of an expanding non-central shell- 
crossing singularity. This singularity reaches R = +oo at t = —0, then disappears 
and leaves a regular spacetime. 

3.1.2. Self-similarity of the zeroth kind 

In the case of self-similarity of the zeroth kind, the Einstein equations imply 
that the quantities m and fi must be of the forms 

2Gm = r[Mi(0 + r 2 M 2 (£)], (3-42) 
SttG/x = ^[Wi(0 + r 2 W 2 (0], (3-43) 

where £ = re~ t . In this case, the basic equations are the following: 

Mi + M[ = IFi5 2 (5 + S'), (3-44) 
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Fig. 2. Same as Fig. 1, but for 1 < a < 3/2. Here, r = r-2 corresponds to the extremum of AH2. 
The trajectory of AH2 is ingoing (outgoing) null at r = n in the left (right) figure. 



3M 2 + M 2 = W 2 S 2 (S + S'), (3-45) 

M{ = 0, (3-46) 

M' 2 = 0, (3-47) 

Mi = S[l - e~ 2 * {S + S')\ (3-48) 

M 2 = 55 /2 Cq 2 , (3-49) 

W[S = -Wi{V'S + 2S'), (3-50) 

WzS = -W 2 (&'S + 2S'), (3-51) 

S" + 5' = (5 + S')&', (3-52) 

S'(S' + 2W'S) = W 2 S 2 c 2 , (3-53) 
25(5" + 25') - 2V'S(S + 5') = -S' 2 - S 2 + e 2 *(l - Wi5 2 ), (3-54) 

255" + 5' 2 = 0, (3-55) 

(5 + 5') 2 = e 2 ^, (3-56) 



where the prime denotes the differentiation with respect to ln£. Equation (3-56) 
gives exp(2<P r ) = (5 + 5') 2 , so that Mi = is concluded from Eq. (3-48). The 
equality M\ = implies W\ = 0, from Eq. (3-44), because |5 + 5'| = exp(lP r ) cannot 
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Region I (Rr>0, U >0) 

Region II (Rr<0,/i <0:unphysical) 

Region III (Rr>0, U >0) 




AH1 

Fig. 3. Same as Fig. 1, but for < a < 1. Here, r = r2 corresponds to the extremum of AH2. The 
trajectory of AH2 is ingoing (outgoing) null at r — n in the left (right) figure. 



be zero. Equation (3-55) gives the evolution equation for S, 

2y' + 3y 2 = 0, (3-57) 

where y = S' / S. This equation can be integrated to yield 

S = |si -soln^l 273 , (3-58) 

where sq and si are integration constants. The relation sq ^ must be satisfied for 
the solution to be a non- vacuum solution. The equality Mi = SS' 2 jc\ is obtained 
from Eq. (3-49), while W 2 = 2,y 2 /[cl{l + y)] is obtained from Eqs. (3-53), (3-52) and 
(3-55). The resulting solution is 

e 2 * = c 2 Q , (3-59) 
2 # (si - 2s /3 - s ln£) 2 



e 



( Sl -s lnO 2 / 3 
(si-2s /3-so(lnr-t)) 2 



(3-60) 



|si — so(ln r — t)| 2 / 3 
R = r\s 1 -s ln^| 2/3 

= r|si-s (lnr-t)| 2/3 , (3-61) 
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Region I (Rr>0, U >0) 

Region II (Rr<0, U <0:unphysical) 



] Region III (Rr>0,//>0) 




Fig. 4. Same as Fig. 1, but for a < 0. r = r\ corresponds to the extremum of AH1. The trajectory 
of AH1 is outgoing (ingoing) null at r = r% in the left (right) figure. 

4s 2 



3cq (si - s In £) (si - 2s /3 - s In £) 



4s 2 

4S ° (3-62) 



3cq(si - s (lnr - t))(si - 2s /3 - s (lnr - t))' 
4sr,r 3 



2Gm = ^ r . (3-63) 
9c 2 



With the coordinate transformations t — > i = cot + cosi/so + (2cq/3) ln(|so|/co) and 
r — > f = (|so|/co) 2//3 r, a one-parameter family of solutions parametrized by en, 

ds 2 = .dt 2 + (t-2co/3-c lnf) 2 rf _ 2 + _ 2 Colnf |4/3 d/? 2 ) (3 . 64) 
|i — Co lnrp/ d 
4 

87tG » = 3( t - co In r)(i - 2c /3 - c In f) ' (3 " 65) 

2Gm = ^f 3 , (3-66) 

is obtained. This solution was obtained previously by Benoit and Coley. 35 ) (See also 
Ref. 37).) This solution belongs to the Lemaitre-Tolman-Bondi family of solutions for 
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the marginally bound case. We call this solution the zeroth-kind KSS-LTB solution. 
We omit the bars on i and f hereafter for simplicity. For this solution, the energy 
density is positive in the regions satisfying 

t<c lnr, c lnr + 2c /3 < t. (3-67) 

The singular surfaces are t = colnr and t = colnr + 2co/3. The former is a shell- 
focusing singularity, because the circumferential radius R is zero there, while the 
latter is a shell-crossing singularity, because the quantity 



dR_i {-t + c lnr + (2c /3)}/(-t + colnr) 1 / 3 for t < c lnr, 
dr ~ 1 {t- c lnr- (2c /3)}/(i - c lnr) 1 / 3 for i>c lnr, 



(3-68) 



is zero there. The fluid is collapsing in the regions in which t < Cq In r and expanding 
in the regions in which t > colnr, because the quantity 

OR 2 

-~r = sign(t - c lnr)-r|t - c lnrp 1/3 (3-69) 

is negative in the former case and positive in the latter case. The trapped region is 
given by 

c lnr - (8/27)r 3 < t < c \nr (3-70) 

and 

c lnr < t < c lnr + (8/27)r 3 . (3-71) 

The apparent horizons are 

t = c lnr + (8/27)r 3 , (AH1) (3-72) 

and 

t = c lnr- (8/27)r 3 . (AH2) (3-73) 

Here we consider the trajectories of the apparent horizons. We consider only the 
solution with \i > in which R r > are satisfied. From Eqs. (3-72) and (3-73), we 
obtain 

dt_ j Co /r + (8/9)r 2 for AH1, 

dr AH \ c /r - (8/9)r 2 for AH2. 1 > 

AH2 has its extremum at r = r^ = (9C0/8) 1 / 3 , while AH1 does not. For the null 
geodesies on the apparent horizon, the relation 



dt 
dr 



\ ±|-c /r+(4/9)r 2 | for AH 1, 
± "I ±| -c Q /r- (4/9)r 2 | for AH2 [6 ' <0) 



are satisfied. For r > 0, the trajectory of AH1 is timelike, while the trajectory of 
AH2 is spacelike (timelike) for 

r 3 < (>)(9/2)c . (3-76) 
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The trajectory of AH2 is ingoing null for 

r 3 = (9/2)cq = 4. 



(3-77) 



Next, we consider r = 0. The inequality [i < is satisfied for r — > along AHl, 
which is unphysical, and therefore we consider only AH2. For AH2, the relation 



fdt 




/dt 




\dr 


ah/ 


dr 





±1 



(3-78) 



is satisfied, so that AH2 is outgoing null at r = 0. 



Region I (Rr>0, U >0) 

Region II (Rr<0, U <0:unphysical) 

Region III (Rr>0, U X>) 
t 




SFS 



AH2 



Fig. 5. The figure on the left is for the zeroth-kind KSS-LTB solution [Eqs (3-64)-(3-66)], and that 
on the right is for its time-reversal solution. AH, SFS and SCS indicate an apparent horizon, a 
shell-focusing singularity and a shell-crossing singularity, respectively. Each diagram is divided 
by SFS and SCS into three regions, which we call regions I, II and III. Region II is unphysical, 
because /j, < there. A solid curve represents an apparent horizon whose trajectory is timelike, 
while a dashed curve represents an apparent horizon whose trajectory is spacelike, r — rn 
corresponds to the extremum of AH2. The trajectory of AH2 is ingoing (outgoing) null at 
r — r§ in the left (right) figure. 

A schematic figure for this solution is given in Fig. 5, along with that for its 
time-reversal solution. The solution is unphysical in region II, in which the energy 
density is negative. Region I on the left in Fig. 5 represents a growing black hole, 
while that in its time-reversal solution represents a decaying white hole. Region 
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III on the left in Fig. 5 represents the inner region of an expanding non-central 
shell-crossing singularity that starts from R = +0 at t = — oo, while that in its 
time-reversal solution represents the inner region of a contracting non-central shell- 
crossing singularity that starts from R = +oo at t = — oo. 

3.1.3. Self-similarity of the infinite kind 

In the case of self-similarity of the infinite kind, the Einstein equations imply 
that the quantities m and /x must be of the forms 

2Gm = Mi(0/t 2 + M 2 (0, (3-79) 
8ttG/x = Wi(0/t 2 + W 2 (0, (3-80) 

where £ = r/t. In this case, the basic equations are the following: 



M[ 


= WiS^S', 


(3-81) 


M' 2 


= W 2 S 2 S'. 


(3-82) 


2Mi + M[ 


= 0, 


(3-83) 


M' 2 


= 0, 


(3-84) 


Mi 


— Cq S S 5 


(3-85) 


M 2 


= 5(1 - e- 2 ^' 2 ), 


(3-86) 


(2Wi + Wj)5 


= -W l (&'S + 2S'), 


(3-87) 




= -W 2 (&'S + 2S'), 


(3-88) 


S" 


= SV, 


(3-89) 




= 5'(2lP- / 5 + 5 / ), 


(3-90) 


1 _ w^V* 


= 2SS" + 5 /2 - 2if / 5 / 5, 


(3-91) 


e 2 * 


= S' 2 , 


(3-92) 





= 2S"S + 2SS' + S' 2 , 


(3-93) 



where the prime denotes the differentiation with respect to ln£. Here we have im- 
plicitly assumed t > 0. When we define £ = r/(—t) for t < 0, equations identical to 
Eqs. (3-81)-(3-93) are obtained. We extend solutions analytically to the region in 
which t < if they are analytic at t = 0. Such extended solutions can be obtained 
when we allow solutions to be defined for t < 0. 

Equation (3-92) gives exp(2'P r ) = S' 2 , so that M 2 = is concluded from Eq. 
(3-86). The equality M 2 = implies W 2 = 0, from Eq. (3-82), because \S'\ = exp(tf) 
cannot be zero. Equation (3-93) gives the evolution equation for S, 

2y > + 3y 2 + 2y = 0, (3-94) 

where y = £'/ S 1 . This equation can be integrated to yield 

S=\s 1 -s r 1 \ 2/3 , (3-95) 

where sq and s\ are integration constants. Here, sq cannot be zero, because this gives 
the contradictory relation exp(2^) = 0. The relation M\ = SS' 2 jc\ is obtained from 
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Eq. (3-85), while W l = -lyjcl is obtained from Eqs. (3-90), (3-89) and (3-93). The 
resulting solution is 

e 2<p = cl (3-96) 

2^/2 ^ 2 -2/~2| a-Ii-2/3 

' g s o r ^ |«i — s C I ' 

= ^s\^-sot\-^, (3-97) 
R = S=\s 1 -s C 1 \ 2/3 

= r -2/3| Sir _ Sot |2/ 3) (3 . 98) 



8ttG// = 



2Gm 



-4g 
3c^t(s\r — sot) ' 

9 C 2 t 2 ' 



(3-99) 



9cgr 2 ' 



(3-100) 



When the time-reversal solution is also considered, we can choose so > and si > 
without loss of generality. With the coordinate transformations t — > i = cot and 
r — > f = (so/co) 2 / 3 r~ 2//3 , a one-parameter family of solutions parametrized by s±, 

d.,* = -dt 2 + P\ Sl f- 3/2 - i\- 2 > : 'df 2 + f 2 \ Sl f- 3 ' 2 - I\ 4 ' 3 dS2 2 , (3401) 

2Gm = ^f 3 , (3-103) 

is obtained. We omit the bars on t and f hereafter for simplicity. This solution is an 
asymptotically flat FRW solution for r — > oo and belongs to the Lemaitre-Tolman- 
Bondi family of solutions for the marginally bound case. We call this solution the 
infinite-kind KSS-LTB solution. It is noted that for si ^ 0, this solution has no 
physical center, and r = corresponds to the R = |si| 2 / 3 . For s± = 0, this solution 
gives the flat FRW solution. For this solution, the energy density is positive in the 
regions satisfying 

t < 0, sir~ 3/2 < t, for si > (3-104) 

and 

t < sir" 3/2 , < t, for si < 0. (3-105) 

The singular surfaces are t = sir~ 3 / 2 and t = 0, which corresponds to R = |si| 2//3 . 
The former is a shell-focusing singularity, because the circumferential radius R is 
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zero there, while the latter is a shell-crossing singularity, because the quantity 

dR _ j (_i)/( air -3/2 _ t) l/3 for t < Sir -3/ 2> 

9r ~ \ */(t - Sl r- 3 / 2 )V3 for t > Sl r^ 2 [6 ' LUb) 

is zero there. The fluid is collapsing in the region in which t < sir" 3 / 2 and expanding 
in the region in which t > sir -3 / 2 , because the quantity 

^ = -signer" 3 / 2 - t)\r\ Sl r-^ 2 - t\'^ (3-107) 

Ob O 

is negative in the former case and positive in the latter case. The trapped region is 
defined by 

Sir -3/2 _ ( 8 /27)r 3 < t < Sl r^/ 2 (3-108) 

and 

sir" 3 / 2 < t < sir" 3 / 2 + (8/27)r 3 . (3-109) 
The apparent horizons are 

t = Sl r~ 3/2 + (8/27)r 3 , (AH1) (3-110) 

and 

t = sir" 3 / 2 - (8/27)r 3 . (AH2) (3-111) 

Here we consider the trajectories of the apparent horizons. We consider only the 
solution with fi > in which R r > are satisfied. From Eqs. (3-110) and (3-111), 
we obtain 



dt 
dr 



J -(3 Sl /2)r- 5 /2 + (8/9)r 2 forA Hl, 
AH I -(3si/2)r- 5 / 2 - (8/9)r 2 for AH2. 1 ' 



AH1 has its extremum at r = rj = (27so/16) 2 / 9 , while AH2 does not. For the null 
geodesies on the apparent horizon, the relation 



dt 
dr 



±\(3s 1 /2)r~ 5 / 2 + (4/9)r 2 | for AH 1 



± | ±|(3si/2)r- 5 / 2 - (4/9)r 2 | for AH2 ^ 3 ' 113 ^ 
are satisfied. In the region r > 0, the trajectory of AH1 is spacelike (timelike) for 

r 9 / 2 < (>)(27/4)si. (3-114) 

It is outgoing null for 

r 9 / 2 = (27/4)5i = rg /2 . (3-115) 
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Region I (Rr>0, fl >0) 

Region II (Rr<0,# <0:unphysical) 

Region III (Rr>0, U >0) 




Fig. 6. The figure on the left is for the infinite-kind KSS-LTB solution [Eqs (3-101)-(3-103)], and 
that on the right is for its time-reversal solution. AH, SFS and SCS indicate an apparent 
horizon, a shell-focusing singularity and a shell-crossing singularity, respectively. Each diagram 
is divided by SFS and SCS into three regions, which we call regions I, II and III. Region II is 
unphysical, because [i < there. A solid curve represents an apparent horizon whose trajectory 
is timelike, while a dashed curve represents an apparent horizon whose trajectory is spacelike. 
r = T7 corresponds to the extremum of AHl. The trajectory of AHl is outgoing (ingoing) null at 
r = r$ in the left (right) figure. It is noted that r — > +0 and r —* +oo correspond to R — » |si| 2//3 
(constant) and R — > +oo, respectively. 



fdt 




/dt 




\dr 


ah/ 


dr 





The trajectory of AH2 is timelike. Next, we consider r = 0. The inequality /i < 
is satisfied for r — > along AH2, which is unphysical. Therefore we consider only 
AHl. For AHl, the relation 

= =Fl (3-116) 

is satisfied, so that AHl is ingoing null at r = 0. 

A schematic figure for this solution is given in Fig. 6, along with that for its 
time-reversal solution. The solution is unphysical in region II, in which the energy 
density is negative. Region I on the left in Fig. 6 represents a decaying white hole, 
while that in its time-reversal solution represents a growing black hole. Region III 
on the left in Fig. 6 represents the formation of a shell-crossing singularity at t = 
that corresponds to R = |si| 2 / 3 from a regular spacetime with no physical center. 
That in its time-reversal solution represents a shell-crossing singularity at t = that 
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corresponds to R = |si| 2//3 and a regular spacetime with no physical center for t > 0. 
3.2. ^ parallel to U» 

In the parallel case for dust solutions of the second, zeroth and infinite kinds, 
<P' = is obtained from basic equations. Therefore we set exp(^) = cq, where cq is 
a positive constant, in this subsection. 

3.2.1. Self-similarity of the second and zeroth kinds 

In these cases, the Einstein equations imply that the quantities m and fi must 
be of the forms 



2Gm = tMi (r) 4- t 3 ~ 2a Mo(r) 


(3- 1171 


o7r(_T/i — i i ^rj -r I v / i / 2 v J • 


1 15 ) 


The basic equations are the following: 




Mi = 0, 


(3-119) 


(3 - 2a)M 2 = 0, 


(3-120) 


M[ = W!S 2 S', 


(3-121) 


M' 2 = W 2 S 2 S', 


(3-122) 


Mi = S[l - S' 2 }, 


(3-123) 


M 2 = Cq 2 S 3 , 


(3-124) 


= -W u 


(3-125) 


= {2a - 3)W 2 , 


(3-126) 


25"5 + 5' 2 = 1 - VFi5 2 , 


(3-127) 


3 = c 2 VF 2 , 


(3-128) 


S' 2 = 1, 


(3-129) 


2« - 3 = 0, 


(3-130) 



where the prime denotes the differentiation with respect to r. Equation (3-130) shows 
that a = 3/2 is the only case, and there are no solutions in the case of self-similarity 
of the zeroth kind. Then, Mi = 0, W± = and W 2 = 3/c 2 , are obtained from Eqs. 
(3-119), (3-125) and (3-128). Equation (3-129) gives the evolution equation for S, 



,/2 



(3-131) 



which is integrated to yield 



S = ±r + s , (3-132) 



where sq is an integration constant. The relation M 2 = S 3 /cq is obtained from Eq. 
(3-124). S must be positive for the positive mass solution. The resulting solution is 
given by 



e 2<p = c 2 , (3-133) 
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R = t(s ±r), (3-134) 

8nG ft = (3-135) 

2Gm = c 2 (s ±r) 3 . (3-136) 

This solution is the flat FRW solution. With the coordinate transformations t — > 
t = (2co/3)t 3 / 2 and r — > f = (3/2co) 2 / 3 (so ± r), the solution in the more usual form 

ds 2 = -dt 2 + t 4/3 (df 2 + r 2 dQ 2 ), (3-137) 
8ttG/x = ^, (3-138) 

2Gm = If 3 (3-139) 

is obtained. 

3.2.2. Self-similarity of the infinite kind 

In this case, the Einstein equations imply that the quantities m and ii must be 
of the forms 

2Gm = M(r), (3-140) 
8ttGh = W(r), (3-141) 

and the Einstein equations and the equations of motion for the matter field (2-4)- 
(2-11) are written 

M = S(l — S' 2 ), (3-142) 
M' = WS'S 2 , (3-143) 
2SS" + S' 2 = 1 - WS 2 , (3-144) 
S' 2 = 1, (3-145) 

where the prime denotes the differentiation with respect to r. Because Eqs. (3-144) 
and (3-145) imply W = 0, there are no solutions in this case. 

3.3. f orthogonal to 

For the dust case, in the case of self-similarity of the second kind, Eqs. (2-4) 
and (2-31) are contradictory, and in the case of self-similarity of the infinite kind, 
Eqs. (2-4) and (2-32) are contradictory, and therefore there are no solutions in these 
cases. 

3.3.1. Self-similarity of the zeroth kind 

In this case, the Einstein equations imply that the quantities m and fx must be 
of the forms 

2Gm = rMi(i) + r 3 M 2 (t), (3-146) 
8irGn = r- 2 W 1 (t) + W 2 {t). (3-147) 
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In this case, the basic equations are the following: 



Mi = 


W X S 3 , 


(3-148) 


3M 2 = 




(3-149) 


Mj = 
i 


0, 


(3-150) 


Mo = 


o 


(3-151) 


Mi = 


5(1 _ e~ m S 2 ), 


(3-152) 


Mo = 


SS' 2 , 


(3-153) 


W[5 = 


-W 1 (VS + 2S'), 


(3-154) 


W^S = 


-W 2 (&'S + 2S'), 


(3-155) 


s' = 


sw, 


(3-156) 


-WxS 2 = 


e~ 2 *S 2 - 1, 


(3-157) 


w 2 s 2 = 




(3-158) 


= 


1 - SV 2 ^, 


(3-159) 


= 


25"5 + 5' 2 , 


(3-160) 



where the prime denotes the differentiation with respect to t. We obtain exp(2^) = 
S 2 from Eq. (3-159), so that M\ = is concluded from Eq. (3-152). The relation 
Mi = implies W\ = 0, from Eq. (3-148). Equation (3-160) gives the evolution 
equation for S, 



2y> + 3y 2 = 0, 

where y = S'/S. This equation can be integrated to yield 

S = \ Sl -s t\ 2 /\ 



(3-161) 



(3-162) 

where so and si are integration constants. The equalities M 2 = SS' 2 and W 2 = 3y 2 
are obtained from Eqs. (3-153) and (3-158), respectively. The resulting solution is 



R 



\si-s t\ 4 ^, 
r\ Sl -s t\ 2 ^, 

2 
5 



4s 2 



3(si - s t) 2 ' 



2Gm = -sir 3 . 



(3-163) 
(3-164) 

(3-165) 
(3-166) 



This solution is the flat FRW solution. With the coordinate transformations t - 
i = si/so — t and r — > f = | So 1 2 / 3 r , the solution in the more usual form 

ds 2 = -dt 2 + i A / 3 {dr 2 + fW), (3-167) 

' (3-168) 



8ttGm= 

4 Q 

2Gm = -f 3 



(3-169) 



is obtained. 
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§4. Corrections and complements to I 

In this section, we give some corrections and complements to I. 39 -* The first is 
only a typographical error: The right-hand side of Eq. (4.44) in I should be replaced 
by a-H'^-^^r' 3 . The others are more important: It is found that the solutions 
that we called the dynamical solutions (A) and (B) result in the singular static 
solution and the flat FRW solution, respectively Although the analytic forms of the 
general solutions have not been obtained for the infinite kind in the tilted case with 
p + fi = (Eqs. (3.126)-(3.132) of I), they are found to be the Nariai solutions, 45 ) so 
that the solution, which is a special solution and we called the yl-cylinder solution, is 
also the Nariai solution. We give a detailed demonstrations of them in the following. 

In the case P2 = W2 = and a = 3/2 considered in §3.1.2 of I, Eqs. (3.15) and 
(3.55) there give 

= y[(l + K)(9K + 2)y' + 2(3K + 1)(1 + K) 2 yy' 

+3{3K + 1)(1 + Kfy 2 + 3(1 + K)(9K + 2)y + 3{K 2 + 6K + 1)]. (4-1) 

Eliminating y' from Eq. (4-1) by using Eq. (3.51) of I, we obtain an algebraic equation 
for y: 

= Ky\l + K)y + 1][6(1 + K){3K + l)y + (11 + 15K)}. (4-2) 

Then, because K 7^ has been assumed, y must be constant. From Eq. (3.51), 
the constant y must be equal to —{K — 3)/[3(l — K)(l + K)] or 0. Substituting 
y = —(K - 3)/[3(l - K)(l + K)\ into Eq. (4-2), we obtain 

= K 2 (K-3) 2 (3K + 5)(3i^-l), (4-3) 

and therefore K = 1/3 is concluded. When K = 1/3, which implies y = — 1, 
Eq. (3.55) of I gives the contradictory relation exp(2!^) = 0, and hence y must be 
zero. The equality y = results in the singular static solution given by Eqs. (3.30)- 
(3.35) of I. 

In the case P2 = W2 = and K / — 1 considered in §3.3.2 of I, the relations 

2^ = -^ + S -- S -, (4-4) 
Wi S' s ' y 1 



2$' 



2K 

:(S' 2 + SS") +2(1 + 3K)S'S" 



l + K 



2K 

: SS' + (1 + 3K)S' 2 



l + K 



(4-5) 



are obtained from Eqs. (3.140) and (3.141) of I. Substituting Eqs. (4-4) and (4-5) 
into Eq. (3.115) of I and eliminating y' by using Eq. (3.139) of I, we obtain 

K[3{1 + 3K){1 + K)y 2 + (9K + 5)y + 2] = 0. (4-6) 

Because we have assumed the solution of Eq. (4-6) is y = — 2/[3(l + K)], 

which is consistent with Eq. (3.139) of I for K 7^ 1, while it is consistent with 
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Eq. (3.141) of I for K = 1. This solution is found to be the flat FRW solution: 

e* = co, (4-7) 

e*= - 2 - s C 2mi+K) K (4-8) 
3(1 + K) s v ; 

5 = sor 2/[3(i+*0] ) (4 . 9) 

M 2 = 0, (4-11) 
4K 

Pl = Wl = 3(iT¥F Co " 2 ' (4 - 12) 

P 2 = KW 2 = 0, (4-13) 



where cq and so are positive constants. With the coordinate transformations t — > 
i = cot and r — > f = (so/cj^^^^r -2 /^ 1 " 1 "^, a more usual form of the solution, 

ds 2 = + £4/13(1+*)] (df 2 + f 2^ 2)) (4 . M) 

2Gm= 9(1 + ^^1+*) ' ^ 
4K 

8nGp = SttGK^ = 3(1 + K)2 ^ 2 , (4-16) 

is obtained. 

We show that a 4-dimensional spherically symmetric spacetime, which contains 
a perfect fluid with p + fi = and has a constant circumferential radius, is described 
by the Nariai solution. In the following, Greek indices denote the 4-dimensional 
spacetime, uppercase Latin indices the 2-dimensional (1 + 1) reduced spacetime, and 
lowercase Latin indices the 2-spheres. We write the general spherically symmetric 
spacetime as a manifold M = M 2 x S 2 with the metric 

= diag(#AB, R 2 lab), (4-17) 

where g A B is an arbitrary Lorentzian metric on M 2 , R is a scalar on M 2 and j a b is 
the unit curvature metric on S 2 . With the shorthand 

v A = (4-18) 

where | denotes the covariant derivative on M 2 , the Einstein equations for a perfect 
fluid with p + fx = in spherical symmetry, in the 2 + 2 split, are 

-2(v A \ B + v A v B ) + {2v c lC + 3v c v c - R- 2 )g A B = SttG^ab, (4-19) 

v c lC + v c v c -K = -8vrG^, (4-20) 

where 1Z = (1/2)R A A is the Gauss curvature of g A B- The energy-momentum con- 
servation equations give that fi is constant (see Eqs. (A-2) and (A-3)). When R is 
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constant, Eqs. (4-19) and (4-20) give that R 2 = 1/(8tvG/j,) and U = 8itGfi, i.e., \i 
must be positive and M 2 is a 2-dimensional positive constant curvature spacetime. 
Since a 2-dimensional positive constant curvature spacetime is the 2-dimensional 
de-Sitter spacetime (dS 2 ), the resulting 4-dimensional spacetime is the Nariai space- 
time, which is dS 2 x S 2 with a constant circumferential radius. Since the spherically 
symmetric kinematic self-similar perfect-fluid solutions of the infinite kind in the 
tilted case with p + /j, = 0, which satisfies Eqs. (3.127)-(3.132) of I, have a constant 
circumferential radius, they are found to be the Nariai solutions, 45 ^ although the 
analytic forms of the general solutions in the self-similarity coordinates have not 
been obtained. 

§5. Spherically symmetric kinematic self-similar perfect-fluid solutions 

of the first kind 

5.1. Tilted case 

In this case, complete classification has already been carried out. (See Refs. 26) 
-29).) 

5.2. f parallel to 

In this case, the Einstein equations imply that the quantities m, fi and p must 
be of the forms 

2Gm = tM(r), (5-1) 
8ttG// = t~ 2 W(r), (5-2) 
MGp = t~ 2 P{r). (5-3) 

The resulting equations for a perfect fluid, (2-4)-(2Tl), reduce to the following: 



p = 


>• 


(5-4) 


(P + W)& = 


-P', 


(5-5) 


M = 


-PS 3 , 


(5-6) 


M' = 


WS 2 S', 


(5-7) 


= 


&s, 


(5-8) 


M = 


S(l + e- 2 *S 2 -S >2 ), 


(5-9) 


1 - ws 2 = 


-3S 2 e- 2<p + 2SS" + S' 2 , 


(5-10) 


-1 - PS 2 = 


sV 2# - 2$'S'S - S' 2 , 


(5-11) 



where the prime denotes the differentiation with respect to r. These equations result 
in 

= co, (5-12) 
P = -\w = pq, (5-13) 
M = -p S\ (5-14) 



A Classification of Kinematic Self-Similar Perfect-Fluid Solutions II 27 

S' 2 = 1 + (cq 2 + P0 )S 2 , (5-15) 
S" = (c 2 + Po )S, (5-16) 

where cq is a positive constant, and the constant po must be negative for positive 
energy density. The resulting solution is the FRW solution: 

(1/a) sin (ar ), forco 2 +p <0, (closed FRW ) 

si ± r, for c 2 + p = 0, (flat FRW) (5-17) 
(1/a) sinh (ar), for c^ 2 + p > 0, (open FRW ) 

2Gm = -p S 3 t, SnGp = -8ttG(h/3) = (5-18) 

where s\ is an integration constant, and a = \J —{cq 2 + Po) for the closed FRW 



solution, while a = y c o +Po f° r the open FRW solution. S must be positive for 
a positive mass solution. For the closed and open FRW solutions, the coordinate 
transformations t — > i = cot and r — ► f = ar are used to write the solution in the 
more usual form 

ds 2 = -dt 2 + -^{df 2 + S 2 dQ 2 ), (5-19) 

2Gm = -^rS 3 i, (5-20) 
c a A 

2 

8ttG P = -8vrG(/V3) = ^p, (5-21) 

| sinhf. (the open FRW ) 

For the flat FRW solution, the coordinate transformations t t = cot and r — ► f = 
(r + s\)/cq are used to write the solution in the more usual form 

ds 2 = -dt 2 + i 2 (df 2 + f 2 dQ 2 ), (5-23) 
2Gm = if 3 , (5-24) 

8nGp = -8irG(n/3) = (5-25) 

t z 

In the vacuum case (i.e., that in which po = 0), the Minkowski solution in the 
Milne form, 

ds 2 = -dt 2 + P(df 2 + sinh 2 fdQ 2 ), (5-26) 

is obtained. 

5.3. £ M orthogonal to 

In this case, the Einstein equations imply that the quantities m, \x and p must 
be of the forms 

2Gm = rM{t), (5-27) 
8itGh = r- 2 W(t), (5-28) 
8-irGp = r- 2 P(t). (5-29) 
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It is noted that in this case, the solution is always singular at the physical center. 
The resulting equations for a perfect fluid, (2-4)-(2-ll), reduce to the following: 



P = w, 

-W'S = (P + W)(V'S + 2S'), 
M = WS 3 , 
M' = -PS 2 S', 

= V'S, 
M = S(l + S' 2 -e- 2 *S 2 ), 
1 - WS 2 = -S^ + e'^S 2 , 
-1-PS 2 = 2S"S + S' 2 - Ze- 2 *S 2 , 


(5-30) 
(5-31) 
(5-32) 
(5-33) 
(5-34) 
(5-35) 
(5-36) 
(5-37) 


where the prime denotes the differentiation with respect to t. These equations result 
in 


e = ci, 
M = WS 3 , 
P = W = w S- 4 , 


(5-38) 
(5-39) 
(5-40) 




(5-41) 


where c\ and wq are positive constants. A positive value of wo is required for a 
positive energy density. c\ can be set to unity with the transformations c\r — > r, 
t/ci — > t, 5/ci — > 5 and wq/c 2 — ► tt)o- In this case, a general solution parametrized 
by t«o can be obtained. The function 5 is given by 


S" 2 = l/2±yJwo- l/4sinh(2t + t ) 


(5-42) 


for wo > 1/4, 




S 2 = 1/2 ± sJi/A-wq cosh(2t + t ), or (1 ± VI - 


4w )/2 (5-43) 


for < w < 1/4, and 




S 2 = 1/2 + e ±2t+t \ 1/2 - e ±2t+t0 , or 1/2 


(5-44) 



for wq = 1/4, where to is a constant and can be set to zero with the transformation 
±2t + to — > ±2t. We call this solution the "singular stiff-fluid solution." The inequal- 
ity wq > must be satisfied for a positive energy density. Also, the sign of S must 
be positive for a positive mass solution. The function S in various cases is plotted 
in Fig. 7. For < wq < 1/4, there exist static solutions with singular centers. Next, 
we discuss the dynamical solutions. For wq > 1/4, the dynamical solution represents 
a big-bang cosmological model or gravitational collapse with a singular center. For 
< wq < 1/4, this solution represents a bouncing cosmological solution with a sin- 
gular center or a big-bang cosmological model with a singular center that ends in a 




Fig. 7. The function S of the singular stiff- fluid solution for (a) wo > 1/4, (b) wo = 1/4 and (c) 
< wo < 1/4. 
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big-crunch. For wq = 1/4, the solution S 2 = l/2 + exp(±2t) represents an expanding 
or contracting universe with a singular center, while the solution S 2 = 1/2— exp(±2t) 
represents a big-bang cosmological model or gravitational collapse with a singular 
center. The trapping condition 2Gm < R corresponds to wq > S 2 . 
In the vacuum case, these equations result in 

e* = ci, (5-45) 
S' 2 = -l + c^ 2 S 2 , (5-46) 
S" = c^ 2 S, (5-47) 

where c\ is a constant. The solution is S = c\ cosh(cJ" 1 t), which is the Minkowski 
spacetime with the metric 

ds 2 = -r 2 dt 2 + c\dr 2 + c\r 2 cosh 2 (c^H)df2 2 . (5-48) 

This metric can be obtained with the coordinate transformations (t, r) — ► (r, x), with 
r = c\r sinh(c^f 1 t) and x = c±r cosh(c ] ^ 1 t) from ds 2 = —dr 2 + dx 2 + x 2 dQ 2 . 

§6. Summary and discussion 



Table I. Kinematic self-similar solutions for the dust case. Here, || and _L denote the parallel and 
orthogonal cases, respectively. 



Self-similarity 


Solution 


Equation number 


1st (tilted) 


see Ref. 25) 




1st (||) 


none 




1st (_L) 


none 




2nd (tilted) 


2nd-kind KSS-LTB (see also Refs. 32), 35), 37)) 


(3-23)-(3-25) 


2nd (||) 


fiat FRW {a = 3/2) 




2nd (_L) 


none 




zeroth (tilted) 


zeroth-kind KSS-LTB (see also Refs. 35), 37)) 


(3-64)-(3-66) 


zeroth (||) 


none 




zeroth (_L) 


flat FRW 




infinite (tilted) 


infinite-kind KSS-LTB 


(3101)-(3-103) 


infinite (| ) 


none 




infinite (_L) 


none 





We have classified kinematic self-similar spherically symmetric dust solutions 
of the second, zeroth and infinite kinds and perfect-fluid solutions, together with 
vacuum solutions of the first kind, in which the kinematic self-similarity vector is 
either parallel or orthogonal to the fluid flow. Except for the cases in which there are 
no solutions, the governing equations can be integrated to obtain exact solutions. 
These kinematic self-similar dust solutions belong to the Lemaitre-Tolman-Bondi 
family of solutions only for the marginally bound case, while, in contrast, there 
exist some first-kind dust solutions that belong to this family for the non-marginally 
bound case. We give a proof of this in Appendix A. We also presented corrections 
and complements to I. 39 ^ It is found that the solutions which in I we called the 
dynamical solutions (A) and (B) result in the singular static solution and the flat 
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FRW solution, respectively. It is found that the spherically symmetric kinematic self- 
similar perfect-fluid solutions of the infinite kind in the tilted case with p + fi = 0, 
which satisfies Eqs. (3.127)-(3.132) of I and includes the special solution which in 
I we called the yl-cylinder solution, are the Nariai solutions. 45 -* The results of the 
present work, together with the corrected results of I are listed in Tables I-V. 

We first summarize the results of this work for the first-kind cases. In the parallel 
case, the flat, open and closed FRW solutions with p = — f/,/3 are the only solutions. 
In the orthogonal case, the equation of state must be p = fi, and a new exact solution, 
which we call the singular stiff-fluid solution, is the only solution. This solution is 
singular at the physical center. 

Next, we summarize the results for the dust case. There are no dust solutions 
in the parallel case of the second kind (except for a = 3/2), zeroth kind and infinite 
kind and in the orthogonal case of the second kind and infinite kind. The flat FRW 
solution is the second-kind solution in the parallel case for a = 3/2 and also the 
zeroth-kind solution in the orthogonal case. In the tilted case, general solutions, 
which represent a collapsing or an expanding dust fluid, can be obtained in closed 
form in the cases of the second kind (second-kind KSS-LTB solution), zeroth kind 
(zeroth-kind KSS-LTB solution) and infinite kind (infinite-kind KSS-LTB solution). 
The second-kind KSS-LTB solution with a = 3/2 is the flat FRW solution. The 
second-kind KSS-LTB solution with a / 3/2, the zeroth-kind KSS-LTB solution 
and the infinite-kind KSS-LTB solution are characterized by a single parameter, 
and the second-kind KSS-LTB solution for a / 3/2 and the infinite-kind KSS-LTB 
solution include the flat FRW solution as a special case. In all three solutions, if 
they are not the flat FRW solutions, there exist a shell-crossing singularity and a 
shell-focusing singularity. There exists an unphysical region, in which the energy 
density is negative. We considered only physical regions. In the collapsing regions, 
the second-kind KSS-LTB solutions for < a < 1 and 1 < a < 3/2 represent the 
formation of a shell- focusing singularity from a regular spacetime. For 1 < a < 3/2, 
this singularity is covered; i.e., this solution represents the formation and growth of 
a black hole. It is important to clarify whether the singularity is naked or covered for 
< a < 1. The second-kind KSS-LTB solution for a < and 3/2 < a, the zeroth- 
kind KSS-LTB solution and the infinite-kind KSS-LTB solution represent a growing 
black hole in the collapsing region, although they do not represent the formation of 
a black hole. These solutions could be important in modeling the growth of a black 
hole. The second- kind KSS-LTB solution for < a < 1 and 1 < a < 3/2 in the 
expanding region could be important in modeling the evolution of a cosmic void. 
The infinite-kind KSS-LTB solution with s\ / does not have a physical center. 
This solution could be attached to an inner region that is represented by another 
solution. This solution could also be important in modeling dust collapse or the 
evolution of a cosmic void. In the high-density region, a relativistic polytrope gas 
might be approximated by a dust, so that the central region in the generic collapse 
of a gas could converge to one of these kinematic self-similar solutions. 

Adding the results of the present work to those of Refs. 25), 26), 38), 39), the 
classification of spherically symmetric kinematic self-similar solutions with the vac- 
uum, dust or perfect fluid obeying a relativistic polytropic equation of state or that 
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of the form p = Kfi, has been completed. 



Table II. Kinematic self-similar solutions for an equation of state of the form p = K/j, 1 [EOS(I)], 
where K and 7 are constants. It is assumed that K 7^ and 7 / 0,1. Here, || and _L denote the 
parallel and orthogonal cases, respectively. (See I 39 ' for the values of the parameters.) 



Self-similarity 


Solution 


1st (tilted) 


none 


1st (||) 


none 


1st (_L) 


none 


2nd (tilted) 


none 


2nd (||) 


none 


2nd (_L) 


none 


zeroth (tilted) 


none 


zeroth (| ) 


none 


zeroth (_L) 


flat FRW 


infinite (tilted) 


none 


infinite (| ) 


All static solutions with the EOS (I) 


infinite (_L) 


none 



Table III. Kinematic self-similar solutions for an equation of state of the form that p — Kn 1 and 
/i = rrify-n + p/(j — 1) [EOS(II)], where the constant mt and n(t,r) represent the mean baryon 
mass and the baryon number density, respectively, and K and 7 are constants. It is assumed 
that K and 7 / 0, 1. Here, || and _L denote the parallel and orthogonal cases, respectively. 
(See I 39 ' for the values of the parameters.) 



Self-similarity 


Solution 


1st (tilted) 


none 


1st (||) 


none 


1st (_L) 


none 


2nd (tilted) 


none 


2nd (||) 


closed FRW 


2nd (_L) 


none 


zeroth (tilted) 


none 


zeroth (||) 


none 


zeroth (1.) 


flat FRW 


infinite (tilted) 


none 


infinite (||) 


All static solutions with the EOS (II) 


infinite (1.) 


none 
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Table IV. Kinematic self-similar solutions for p — Kfi [EOS(III)], where if is a constant and 
— 1 < K < 1 (K =fc 0). Here, || and _L denote the parallel and orthogonal cases, respectively. 
(See the main text and I 39 ' for the values of the parameters a and K.) It is noted that the 
kinematic self-similar solution of the infinite kind with K = — 1 for a tilted case is the Nariai 
solution 45 - ) , although the analytic forms of the solutions in the self-similarity coordinates have 
not been obtained. 



S6lf- similarity 


Solution 


Equation number 




q PP Refs 9fi~l 1 9) 




1st (II) 


FRW (K = -1/3) 




1st (_L) 


Singular stiff-fluid (K = 1) 


(5-38)-(5-44) 


2nd (tilted) 


flat FRW 
Singular static 


(3.41)-(3.46) in I 


2nd (||) 


flat FRW 




2nd (_L) 


Singular static 


(5.28)-(5.31) in I 


zeroth (tilted) 


de-Sitter (K = -1) 
Singular static 


(3.96)-(3.101) in I 


zeroth (||) 


de-Sitter (K = -1) 




zeroth (_L) 


flat FRW 

de-Sitter (K = -1) 




infinite (tilted) 


flat FRW 
Nariai (K = -1) 


(3.127)-(3.132) in I 


infinite (||) 


All static solutions with the EOS (III) 




infinite (_L) 


none 





Table V. Kinematic self-similar solutions for the vacuum case. Here, || and _L denote the parallel 
and orthogonal cases, respectively. 



Self-similarity 


Solution 


1st (||) 


Minkowski 


1st (_L) 


Minkowski 


2nd (tilted) 


Minkowski (all a) 




Schwarzschild (a — 3/2) 


2nd (||) 


none 


2nd (_L) 


none 


zeroth (tilted) 


Minkowski 


zeroth (||) 


none 


zeroth (_L) 


Minkowski 


infinite (tilted) 


none 


infinite (| ) 


Minkowski 




Schwarzschild 


infinite (_L) 


none 



Appendix A 

Flatness of 3-Surfaces for Kinematic Self-Similar Dust Solutions 

In this appendix, we give the integrability conditions for kinematic self-similarity 
of the Einstein equations and a proof of the flatness of 3-surfaces for kinematic self- 
similar dust solutions. The covariant derivative of can be decomposed into its 
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irreducible parts as 

Ufi;u = 0> + ^Ohpv + UJ^ U - UpU v , (A-l) 

where 9 = g^O^v (with 6^ v = h^h u x u K .\) is the expansion scalar, a^ v = 6^ v 
(l/2>)6hn U is the shear tensor (from which we obtain a 2 = (l/2)cr„ I/ cr 



■pu m u±±c aiicoi uciiovji wmui »»c uuua,iii w — \ / / fjiis ) ) ^^t^ 

h^ K h u ^ x u K; \ is the vorticity tensor (from which we obtain cj 2 = (1/2)0^0;^), and 
?V = Ufi-,i/U u is the acceleration vector. 44 ) The conservation laws become 

A+(Ai+p)0 = O, (A-2) 

(H + = -p jV h v M , (A-3) 

where the dot represents operation with while the Ricci identities with the 
Einstein equations become 

+ \e 2 - ^ + 2{a 2 - u 2 ) + ^(/i + 3p) = 0, (A-4) 

(V + = a^u" + ^"V^a, (A-5) 

(^ ;„ - o>» .„ + +{u? v + o\)u» = 0, (A-6) 

uA M = 2uAv, (A-7) 
fyu A ^A - hp K h u A tt( K;A) - iipU u + u^uv + cr^cr" v 

+ \0apv + V ("^ 2 " \° 2 + ; «) + ^ = °> (A-8) 

fT^ = 2u ( ^ } - ^ A ^ C( W(A "^ + a ( / ; > cw <A (A-9) 



with 



Ep V = Cp KuX u K u\ (A- 10) 

Hp V = \ %K <>°C pavX u«u x , (A-ll) 



where C^ VK \ and f]p VK \ are the Weyl tensor and the totally-skew pseudotensor, re- 
spectively. In the case of zero vorticity (o^j, = 0), the Gauss-Codacci equations 
become 



3 Rfj,u = {-0a K \ - a K \ + U( K -\))hp K h v A 

+u li u v + i V (-^ 2 + 2 ^ 2 + 16 ^ G /" " ^ P ;p) > (A-12) 
3 i? = --0 2 + 2<r 2 + 167rG/i, (A-13) 

where 3 i?^^ is the Ricci tensor of the 3-spaces orthogonal to u^, and 3 i? is the 
corresponding Ricci scalar. 
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Now we derive the integrability conditions for kinematic self-similarity on the 
Einstein equations, i.e., 



= ^GC^ [r pu - , (A- 14) 



where R pv is the Ricci tensor and £ M is a kinematic self-similarity vector that is 
defined by Eqs. (2-18) and (2-19). From Eqs. (2-18) and (2-19), we obtain 

C^g pv = 28g iiV + 2(5 - a)u p u u . (A- 15) 

Hence, we obtain 

L(R pv = (C^r^ u ). K - (C^ K ). V 

= (5 -a) [2a ^ + 29a pv + 2a pK cu K „ + 2a UK uj K ^ + Auj^ k uj k v 

+\og pv (9 + 9) + u p (2co KU < K - 2a VR u n - 2u VK ix R ) 



+u u (2to Kp ' K - 2a pK ii K - 2aj pK u K 
+ Ulx u v {^9(9 + 9)- 2u K -> K 



(A-16) 



where 



W vX = \sr[{Ci9*»);X + {^9nx);u ~ (C^vx);*]- (A-17) 

Decomposing this equation in various ways, we obtain 

(Ct-R pv )u^u v = (5- a)(-8co 2 - 2ii R ;K ), (A- 18) 

(L i R llv )h ilv = 2(8 - a)(9 + 9 2 - 4u 2 ), (A- 19) 
(CtR^uW x = (5- a)(2oj x ^ + 2uj kX ' k 

-4u 2 u x ), (A-20) 

(C^) (V x h u p - -h Xp h^ = 2(5 - a) (a Xp - u p a Xv ii v - u x a m u» 

+9a Xp + a XK uj K p + a pK uj K x 



In the perfect-fluid case, the relation 



+2lu x K u Rp + -h Xp LO z ) . (A-21) 



T pv - -g^T = -(// + 3p)u p u u + - ((j, - p)h pv (A-22) 



is satisfied. Then, we obtain 



1 



u p u u 



+^{£^ + 25fi- £^p-25p}h pu . (A-23) 
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Decomposing this equation, we obtain 



:, ' //(* .y'll'"^ 



; (%i + 2a/i) 



--(%> + 2ap), (A-24) 



h\h\--h Xp h^ 



--(C (P + 25p), (A-25) 
0, (A-26) 
0. (A-27) 



Hence, from Eqs. (A-18)-(A-21) and Eqs. (A-24)-(A-27), the integrability conditions 
for the kinematic self-similarity (A- 14) are obtained as follows: 



(5-a)(-8uj 2 -2u K ' K ) = 8ttG 

2(5 - a)(6 + 6 2 - 4w 2 ) = 8vrG 
2u x ^ + 2u kX ' k -4u; 2 u a = 0, 



-{jC^ + 2a f i) + -(^p + 2ap) 
^(C^ + 2Sfi)-^p + 28p) 



,(A28) 

, (A-29) 
(A-30) 



and 



<?\p ~ u p o\ u u - u x a pp ,w 

+ 0a Xp + o Xk uj k p + a pK u R A + 2uj x K u Kp + Ui Xp u 2 = 0. (A-31) 

These equations must be satisfied in addition to the basic equations. The integrabil- 
ity conditions of the finite kinds (i.e., the first, second and zeroth kinds), assuming 
'physical' self-similarity for the perfect fluid, i.e., 



C^fi = afj,, L(j> = bp, 



(A-32) 



where a and b are constants, were obtained by Coley. 34 ) 

If it is assumed that the fluid is irrotational, so that the vorticity is zero (i.e., 
= 0), comoving coordinates can be adopted, and we obtain 



-6 2 -u». a 
3 ^ 



2a 2 + 



fi + (n + p)0 = O, 
8ttG, 

— (/i + 3p) = 0, 



(A-33) 
(A-34) 

(A-35) 
(A-36) 
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from Eqs. (A-2), (A-3), (A-4), (A-6) and (A-8), respectively, while 
-2(5 - a)u K ;K = 8vrG 
2(5 - a)(9 + 9 2 ) = 8itG 



^C^(i + u(i + ^C^p + 3ap 



3 3 
-C^n + 38(1 - -C^p - 35p 



= <T\ P - u p a\ u u u - u x a Pfl u^ + do\ p 



(A-37) 

(A-38) 
(A-39) 



from Eqs. (A-28), (A-29) and (A-31), respectively. The Gauss-Codacci equations 
must be satisfied in addition. The remaining non-trivial equations serve to define 
and H pu . The Einstien tensor on the hypersurface orthogonal to are obtained 
from Eqs. (A-8), (A- 12), (A- 13) and (A-39) as 



1. 



3 R/j,u — -h^R — E pu + a UK a K 



1 2 
[ik u u ~ -Oa^ - -a h pv . (A-40) 



Now we demonstrate the flatness of 3-surfaces orthogonal to the fluid flow for 
kinematic self-similar dust solutions. When p = 0, Eq. (A-3) gives = 0; i.e., the 
dust fluid is geodesic. From Eq. (A-4), we obtain 



8ttG 



1 



-2a 2 . 



(A-41) 



When (i > 0, 9 cannot be identically zero, since a 2 > 0. Hence the following Lemma 
is shown. 



Lemma 1 If a spacetime that satisfies the Einstein equations contains an irrota- 
tional dust fluid with positive energy density, then 9 cannot be identically zero. 

When p = 0, Eq. (A-37) gives C^(i = —2a(i. Hence, the following Lemma is shown. 

Lemma 2 If a spacetime satisfies the Einstein equations, contains an irrotational 
dust fluid, and admits a kinematic self-similarity vector field of the second, zeroth or 
infinite kind, then the dust fluid exhibits physical self- similarity with a = —2a. 

We now show the following theorem. 

Theorem 1 Any three surface orthogonal to the fluid flow in a kinemetic self-similar 
solution to the Einstein equations of the second, zeroth or infinite kind that contains 
only irrotational dust with positive energy density as a matter field is flat. 

Proof. Equation (A-39) gives 

6-pv + 0<v = 0. (A-42) 
Multiplying this equation by a^ v ', we obtain 

a[& + 9a] = 0. (A-43) 
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Equations (A-37) and (A-35) with Lemma 2 give 

(5 -a)(e + U 2 + 2a 2 + 8ttG| 



0. 



Differentiating this and using Eqs. (A-43) and (A-44), we obtain 



(5 -a) 



o, 



(A-44) 



(A-45) 



which yields, with use of Eq. (A-33) and the equation that is obtained by differenti- 
ating Eq. (A-38) with Lemma 2, 



(5 - a)9 



' : (6 + 6 2 )-8irGfi 



By Lemma 1, this equation reduces to 



(5 -a) 



;(e + 6 2 )-8irGfi 



0. 



= 0. 



Hence we obtain 



2(5 - a) (--e 2 + a 2 + 8irGfjtj = 



(A-46) 



(A-47) 



(A-48) 



from Eqs. (A-44) and (A-47). Finally from Eqs. (A- 12), (A-42) and (A-48), we 
consequently find that 



(5 - afR 



/IV 



0; 



(A-49) 



i.e., 3-spaces orthogonal to are flat in the case of kinematic self-similarity of the 
second, zeroth and infinite kinds, because a / 5. ■ 



If the spacetime is spherically symmetric, the solution can be obtained as the Lemaitre- 
Tolman-Bondi family of solutions. Because the fact that the 3-surfaces orthogonal 
to the fluid flow are flat corresponds to the fact that these solutions are those for 
the marginally bound case, we can show the following corollary. 



Corollary 1 A spherically symmetric kinemetic self-similar solution to the Einstein 
equations of the second, zeroth or infinite kind that contains only dust with positive 
energy density as its matter field is described by the marginally bound Lemaitre- 
Tolman-Bondi family of solutions. 
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